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Abstract 

We define a class of quantum systems called regular quantum graphs. Al- 
though their dynamics is chaotic in the classical limit with positive topological 
entropy, the spectrum of regular quantum graphs is explicitly computable an- 
alytically and exactly, state by state, by means of periodic orbit expansions. 
We prove analytically that the periodic orbit series exist and converge to the 
correct spectral eigenvalues. We investigate the convergence properties of the 
periodic orbit series and prove rigorously that both conditionally convergent 
and absolutely convergent cases can be found. We compare the periodic orbit 
expansion technique with Lagrange's inversion formula. While both meth- 
ods work and yield exact results, the periodic orbit expansion technique has 
conceptual value since all the terms in the expansion have direct physical 
meaning and higher order corrections are obtained according to physically 
obvious rules. In addition our periodic orbit expansions provide explicit ana- 
lytical solutions for many classic text-book examples of quantum mechanics 

that previously could only be solved using graphical or numerical techniques. 
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I. INTRODUCTION 



At a first glance it may seem surprising that many chaotic dynamical systems have 
explicit analytical solutions. But many examples are readily at hand. The shift map |l],|2| 

x n+ i = (2x n ) modi, x n e R, n = 0,1,2,..., (1.1) 

for instance, is "Bernoulli" |1, the strongest form of chaos. Nevertheless the shift map is 
readily solved explicitly, 

x n = (2 n x ) modi, x n G R, n = 0,1,2, (1.2) 

Another example is provided by the logistic mapping 

x n+1 = fj,x n (l - x n ), x n e[0, 1], < n < 4, n = 0,1,2,..., (1.3) 



widely used in population dynamics |T]-§[. For /i = 4 this mapping is equivalent with the 
shift map Q and therefore again completely chaotic. Yet an explicit solution, valid at /i = 4, 
is given by H,§: 



x n = sin 2 (2 n arcsin sfx^) , x E [0, 1]. (1.4) 

Therefore, as far as classical chaos is concerned, there is no basis for the belief that classically 
chaotic systems do not allow for explicit analytical solutions. But what about the quantized 
versions of classically chaotic systems, commonly known as quantum chaotic systems [§,0? 
Here, too, the answer is affirmative. It was shown recently l8Hl2|| that regular quantum 
graphs provide a large class of explicitly solvable quantum chaotic systems. In order to 



strengthen first analytical and numerical results presented in [p|-[T2]|, it is the purpose of 



this paper to show that the explicit periodic orbit expansions obtained in [p|-[T2]| are more 
than formal identities. We will prove below that (i) the spectrum of regular quantum 
graphs is computable explicitly and analytically, state by state, via convergent periodic orbit 
expansions and (ii) the periodic orbit series converge to the correct spectral eigenvalues. 
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The main body of this paper is organized as follows. In Sec. II we summarize the 
basics of quantum graph theory and derive the general form of the spectral equation. In 
Sec. Ill we define regular quantum graphs. In Sec. IV we present the explicit periodic 
orbit expansions of individual eigenvalues of regular quantum graphs. We also specify a 
summation scheme that guarantees convergence to the correct results. The derivations 
presented in Sec. IV are mathematically rigorous except for one step where we interchange 
integration and summation to arrive at our final results. This is a "dangerous operation", 
which is not allowed without further investigation. This point is resolved in Sec. V, where 
we present the analytical proofs that justify the interchange of integration and summation 
performed in Sec. IV. This result establishes that the periodic orbit expansions investigated 
in this paper converge in the usual sense of elementary analysis. In Sec. VI we investigate the 
convergence properties of the periodic orbit series obtained in Sec. IV. We prove analytically 
that there exists at least one quantum graph for which the convergence is not absolute, but 



conditional. According to Riemann's well-known reordering theorem (see, e.g. [I3]|, volume 
II, page 33), it is possible to reorder the terms of a conditionally convergent sum in such 
a way that it converges to any prescribed number. This is the reason why in Sec. IV we 
place so much emphasis on specifying a summation scheme that guarantees convergence 
of the periodic orbit sums to the correct spectral eigenvalues. In Sec. VII we present an 
alternative way of solving the spectral equations of regular quantum graphs: Lagrange's 
inversion formula. Although a perfectly sound and fast-converging method for solving the 
spectral equation, it lacks the physical appeal of periodic orbit expansions which are based 
on concrete physical objects such as periodic orbits and their geometrical and dynamical 
properties. In Sec. VIII we discuss our results, point out promising research directions and 
conclude the paper. 

Since many scientists will find the existence of convergent periodic orbit expansions 
surprising, we found it necessary to establish this result with an extra degree of rigor. This 
necessarily requires strong, but somewhat technical proofs. However, in order not to break 
the flow of the paper, we adopt a hierarchical approach presenting the material in three 



stages. The most pertinent aspects of our results are presented in the main text as outlined 
above. Supporting higher-level, but formal, material is presented in Appendix A. Lower 
level material, such as formulae, definitions and lemmas, is relegated to Appendix B. 

In order for the proofs to be convincing, and to be accessible to a wide readership, we 
used only concepts of elementary undergraduate analysis in our proofs, altogether avoiding 
advanced mathematical concepts, such as distributions fl4]]l5l| . This is an important point 
since our paper seeks mathematical rigor on a common basis acceptable to all readers. In 
this spirit the alert reader will notice that we completely avoid the use of Dirac's delta 



"function" ||16|| . This is necessary since the delta "function" is a distribution fl4| , [i5|1 , a 
concept we found to be highly confusing for the general reader. Although there is nothing 
"wrong" with the delta "function" , if treated properly as a distribution or linear functional, 



the confusion surrounding the delta "function" started with von Neumann's critique [17 



at a time when the modern tools of distribution theory were not yet available. Therefore, 



although ultimately completely unjustified, von Neumann's critique |T7| tainted the delta 
"function" to such an extent that it still can't be used in a rigorous context without causing 
heated debate. Thus, for the sake of simplicity and clarity of our arguments, we prefer to 
avoid it. As a consequence, the presentation of the material in this paper and all of our 
proofs are conducted without using the concept of level densities, which are usually defined 
with the help of Dirac's delta "function". 

We emphasize that there is a fundamental difference between an analytical solution and 
an explicit analytical solution. While the spectral equation for quantum graphs is known 
in great detail ||18|| , and periodic orbit expansions for the spectral density and the staircase 
function of quantum graphs are known |Tj|, these results are all implicit. This means that 
they do not yield the spectral eigenvalues in the form u E n = ...". It is the novelty of 
this paper to obtain explicit analytical expressions of this type for a wide class of quantum 
chaotic systems, and prove their validity and convergence with mathematical rigor. 
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II. QUANTUM GRAPHS AND SPECTRAL EQUATION 



The properties of the Laplacian operator on graphs have been studied in great detail 
in the mathematical literature [1S,2C] and the study of quantum mechanics on graphs has 



attracted considerable attention among chemists and physicists (see, e.g., [21-23] and ref- 
erences therein), especially in the quantum chaos community |T^J2^J2^-^| . The purpose 
of this section is to acquaint the reader with the main ideas of quantum graph theory. 
Since many excellent publications on the theory of quantum graphs are available (see, e.g., 
|T%| , |22| -P7[ ) , we will present only those ideas and facts that are of direct relevance to the 
subject of this paper. 

A quantum graph consists of a network of bonds and vertices with a quantum particle 
travelling on it. An example of a graph with ten bonds and six vertices is shown in Fig. 1. 
The number of bonds is denoted by Nb, the number of vertices by Ny. In this paper we 
focus entirely on finite quantum graphs, i.e. N B , N v < oo. We define directed bonds on the 
graph such that the bond Bij connecting vertex number i with vertex number j is different 
from the bond Bji connecting the vertices in the opposite direction. There are 2Nb directed 
bonds. It is useful to define the linearized bond index A = 1, . . . , Nb'- (ij) \—>A,i<j,i and 
j connected. The index A labels sequentially all Nb directed bonds of the graph with i < j. 
For the directed bonds B^ with i > j we define A(ij) = —A(ji). This way the sign of the 
counting index reflects the directionality of the bond. The network of bonds and vertices 
defines the graph's topology. The topology of a graph alone does not completely specify the 
quantum graph. This is so because, for instance, the bonds of the graph may be dressed 
with potentials 0. Since the quantum graphs we study in this paper are finite, bounded 
systems, their spectra are discrete. 

The spectrum of a quantum graph is obtained by solving the one-dimensional Schrodinger 
equation on the graph subjected to the usual boundary conditions of continuity and quantum 
flux conservation. A particularly useful way of obtaining the spectral equation for a given 
quantum graph is the scattering quantization approach which yields the spectral 
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equation in the form 

det[l - S(k)] — 0, (2.1) 



where S is the quantum scattering matrix of the graph and k is the wave number, 

related to the energy via E = k 2 . For our purposes it is sufficient to know that the S matrix 
is of dimension 2Nb x 2Nb and can be decomposed into [p^ ,|l8|] 



S{k) = D{k)U{k), (2.2) 

where U(k) is a 2Nb x 2Nb unitary matrix and D(k) is a diagonal matrix of the form 

D AA ,(k) =exp[iL A (k)] 5 A)A ,, A,A' = ±l,...,±iV B (2.3) 

with L A (k) e R, L\(k) = L-\(k) > 0. The ordering of the matrices D and U in (|2.2Q is 
neither unique nor important in the present context. It depends on the details of how "in" 
and "out" channels are assigned to the indices of the S matrix. Since the computation of 
the energy spectrum involves only traces and determinants, the precise ordering of D and 
U in ( |2.2| ) does not affect our final results. 

Physically the quantities L\(k) are the time-reversal invariant parts of the bond actions 
|10| , [T2|JT8| . A possible time- reversal breaking part of the bond actions is understood to be 
absorbed in the matrix U |12j| . For simplicity we will in the following refer to L\ as the 
bond action of the bond B A . 

In this paper we focus exclusively on scaling quantum graphs PHT^,^ . In this case the 
matrix U is a constant matrix, independent of k, and the actions L\(k) in (|2.3|) split into 
the product L\(k) = L A k, where G R is a constant, the reduced bond action of the 
bond B A . Physically the scaling case is an important case since it describes systems free 
from phase-space metamorphoses PD[ . Scaling systems of this type arise in many physi- 
cal contexts, for instance in microwave cavities partially filled with a dielectric pl|-|33|, or 
Rydberg atoms in external fields p| . 
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It is possible to write (|2.1| ) as a linear combination of trigonometric functions whose 
frequencies are directly related to the bond actions. In order to derive this representation 
we start by noting that 

2N B 



det[l - S] = det (S 1/2 ) (-4)^ J] sin ( ff //2), 



(2.4) 



1=1 



where <j\ G R are the eigenphases of the unitary matrix S. /,From ( |2.4| ) we obtain the 
important result that 



det(V 1/2 )det[l-S] G R. 



(2.5) 



According to ( p.2|) we have 



det(,S- 1/2 ) = exp(-ir/2) exp I -i E Lf ] k 



(0), 



(2.6) 



A=l 



where r is the eigenphase of the unitary matrix U, i.e. det(U) = exp(ir). Using the 
decomposition ( |2.2|) of the 5* matrix we write the spectral determinant (|2.1| ) in the form 



det 



-U- l ) + D 



(2.7) 



so that we can directly apply the results in Ref. [p5}[ . According to Ref. [|3"5j| , p. 87, the 
determinant ( |2.7| ) is a polynomial in the 2Nb variables exp(iL^ k), A = ±1, . . . , ±Nb, whose 
coefficients are the principal sub-determinants [3l| of —U^ 1 . Using = together with 
the fact that det(?7" 1 ) is the principal sub-determinant of —U^ 1 of order zero [35), we obtain 



that ( 2/7) is of the form 
det 



(_c/-i) + D \= detiu- 1 ) + E E E V.,^) ex P 

n=l ii,...,i„ Q,(n) 



Ae{«l,...,in} 



(2.8) 



where 1 < %\ < %2 < ■ ■ ■ < i n < Ng, is an integer array of length n containing a "1,2- 
pattern", i.e. orf 1 G {1,2}, and A n . ix _ . j n . a (n) are complex coefficients that can be computed 
from the principal sub-determinants of —U~ l . Because of ( |2.5| ), we have 



det (S- 1/2 ) det(U) det (-U' 1 + D) G R. 
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(2.9) 



Defining 



N B 



LO 



E4" 1 



A=l 



and using ( |2.6| ) and fl2.8| ), we obtain ( j2.9| ) in the form 

Nb 

exp[-i(u k + t/2)} + ^ ^ X ^n,...,^;^) ex P {-^/W,-,^™)^ - r 

n=l ii,...,in ct(™) 



(2.10) 



e R, 
(2.11) 



where 



/? _ V 1 ( n ) T (°) 

Pn;ji,...,j„;a(") ~ W ~ 2^ a A ■ 

Ae{ii,...,i n } 



We define the frequencies 



(2.12) 



^n;ji,...,i n ;o( n ) l/^n;ii,...,j n ;a( n ) I - 



(2.13) 



Because of > f° r a ^ A an d the structure of a*™), the largest frequency in ( |2.13|) is 



ujq defined in ( |2.10| ). We now scan the frequencies LO n . ix i n - a (n) defined in Q2.13|) and collect 



the pairwise different ones into a set Q = {u>o, . . . , ojm }, where M + 1 = \Q\ is the number 



of pairwise different frequencies and lum < lom-\ < ■ ■ ■ < u>o- Since the derivation of ( |2.11| ) 
involved only factoring nonzero terms out of the left-hand side of Q2.1D , the zeros of fl2.11p 
and the zeros of ( |2.1p are identical. Taking the real part of the real quantity ( |2.11p shows 
that ( |2.1D can be written in the form 



cos(a;o^ — 7T7o) — = 0, 



(2.14) 



where 



M 



=J2ai COs(uJik - 7T7i 



(2.15) 



i=l 



and di,ji, i = 0, . . . , M are real constants. In general it is difficult to obtain an explicit 



analytical result for the zeros of ( |2.14j ). There is, however, a subset of quantum graphs 
defined in the following section, that allows us to compute an explicit analytical solution of 



s 



III. REGULAR QUANTUM GRAPHS 



A subset of quantum graphs are regular quantum graphs. They fulfil the condition 



a 



I a *l < X > 



(3.1) 



where the constants cij are the coefficients of the trigonometric functions in ( |2.15 ). Although 
regular quantum graphs are a restricted sub-set of all quantum graphs, they are still quantum 
chaotic with positive topological entropy P,[TT]|-|T2"|| . Because of ( |3~T| ) we have < 1 for 

all k, and the zeros k n of ( [2.14 ) are given by: 



k n k n -\- k n ^ ti 1,2,..., 



(3.2) 



where 



k - — 

UJ 



n + n + - + 7o 



(3.3) 



and 



k r . 



LU 



arccos[$(/c n )] 



7T 



(3.4) 



may be interpreted as the average and fluctuating parts of the zeros of ( |2.14j ), respectively. 
Since ( |2.14| ) is the spectral equation of a physics problem, we only need to study the positive 
solutions of ( p.!4|) . Therefore we introduced the constant fi G Z in (|3.3|) which allows us to 
adjust the counting scheme of zeros in such a way that k\ is the first nonnegative solution of 
( |2.14j ). This is merely a matter of convenience and certainly not a restriction of generality. 



Because of (|3.1|) , the boundedness of the trigonometric functions in (|2.15| ) and the properties 
of the arccos function, the fluctuating part of the zeros is bounded. We have 



\u I < u ■— 

I | _ "-max • 



arccos ( a 

2 v ' 



< 



71 

2u ' 



(3.5) 



Therefore, roots of ( |2.14 ) can only be found in the intervals R n = [k n — k max , k n + k v 
We define 
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7T 

k n = — (^ + ^ + 1 + 70), n = l,2, ... (3.6) 

Wo 

and note in passing that 

k n = (k n -i + L)/2. (3.7) 

We also define the open interval I n = (k n -i, k n ) and its closure I n = [k n -i, k n ). For a < 1 we 
have R n C I n . For a — *■ 1 the root intervals i? n grow in size towards J n , but for any a < 1 
the end points fc n _i and k n of J ra are not roots of (|2.14 ). The boundedness of k n also implies 



the existence of two root-free intervals in /„. They are given by = [k n -i, k n — k mSLK ] and 
= [k n + fc max , k n ). Thus, roots cannot be found in the union of these two intervals, the 
root-free zone F n = F^ U F^ C /„. We also have I n = F n U R n . For an illustration of the 
various intervals defined above, and their relation to each other, see Fig. 2. The intervals 
/„ together with their limiting points k n provide a natural organization of the k axis into a 
periodic structure of root cells. 



We now define x = uj$k — 7770, which transforms ([2.14|) into 



M 



cos(x) — $(a;) = 0, $(x) = cos(pjX — rrvi), (3.8) 

8=1 

where pi = Ui/uo and h>i = 7j — Pi7o- Since, as discussed in Sec. II, u is the largest frequency 
in (|2.14|) , we have pi < 1, i = 1, . . . , M, and theorem T2 (Appendix A) is applicable. It 



states that there is exactly one zero x n of (|3.8j ) in every open interval (nir, {n + l)vr), hgZ. 
Consulting Fig. 3 this fact is intuitively clear since the cos function in (|3.8|) is "fast", and 
$(x), containing only frequencies smaller than 1, is a "slow" function. Thus, as illustrated 
in Fig. 3, and proved rigorously by T2 (Appendix A), there is one and only one intersection 
between the fast cos function and the slow $ function in every x interval of length tt. 
Transformed back to the variable k this implies that there is exactly one zero k n of ( |2.14| ) 



in every interval I n . Since this zero cannot be found in the root-free zone F n , it has to be 
located in R n . Thus there is exactly one root k n of ( j2.14|) in every root-interval R n . This fact 
is the key for obtaining explicit analytical solutions of fl2.14| ) as discussed in the following 
section. 
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IV. PERIODIC ORBIT EXPANSIONS FOR INDIVIDUAL SPECTRAL POINTS 



For the zeros of ( p. 14 ) we define the spectral staircase 

oo 

where 



(4.1) 



0, for x < 0, 

0(x) = { 1/2, for x = 0, (4.2) 
VI, for x > 0, x G R 
is Heavyside's function. Based on the scattering quantization approach it was shown 
elsewhere |I8fl that 



1 



N(k) = N(k) + -ImTr ]T - S (k), 



7T 



(4.3) 



i=i 



where 



N(k) = — - (/i+l + 7o), 



7T 



(4.4) 



and S(k) is the unitary scattering matrix ( |2.2| ) of the quantum graph. Since, according to our 
assumptions, S is a finite, unitary matrix, existence and convergence of (|4.3|) is guaranteed 
according to L17, L18 and L19 (Appendix B). Therefore, N(k) is well-defined for all k. 
Since S{k) can easily be constructed for any given quantum graph [|T^,|TS|], ( |4.3| ) provides an 



explicit formula for the staircase function (|4.1| ). Combined with the spectral properties of 
regular quantum graphs discussed in Sec. Ill, this expression now enables us to explicitly 
compute the zeros of fl2.14|) . 



In Sec. Ill we proved that exactly one zero k n of fl2.14p is located in I n = (k n -x,k n ). 
Integrating N(k) from k n -i to k n and taking into account that N(k) jumps by one unit at 
k = k n (see illustration in Fig. 4), we obtain 



N(k) dk = N(k n -i) [k n - k n -i\ + N(k n ) [kn - k n ] . 



Solving for k n and using N(k n -i) = n — 1 and N(k n ) = n (see Fig. 4), we obtain 



(4.5) 
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71 



k n = — (2n + n + 70 ) - /. N(k)dk. (4.6) 

^0 •'fcn 



n — 1 



Since we know N(k) explicitly, (4J5) allows us to compute every zero of ( 2.14 ) explicitly and 



individually for any choice of n. The representation ( |4.6|) requires no further proof since, 
as mentioned above, N(k) is well-defined everywhere, and is Riemann-integrable over any 
finite interval of k. 

Another useful representation of k n is obtained by substituting (|4.3|) with (|4.4j) into (|4.6|) 



and using (|3~3|): 

k n = k n ImTr / Y.l S \ k ) dk - ( 4 - 7 ) 

71 Jk n -i t 

According to theorem T3 (Appendix A) presented in Sec. V, it is possible to interchange 
integration and summation in ( ^.7| ) and we arrive at 

k n = k n ImTr E 7 /- S\k)dk. (4.8) 

In many cases the integral over S l (k) can be performed explicitly, which yields explicit 
representations for k n . 

Finally we discuss explicit representations of k n in terms of periodic orbits. Based on 



the product form ( |2.2| ) of the S matrix and the explicit representation ( |2.3| ) of the matrix 
elements of D, the trace of S(k) 1 is of the form 

Tr S(k) 1 = £ D jldl U juj2 D j2j2 U j2j3 . . . D jhjl U jldl = £ A m [l] exp {iL® I k} , (4.9) 

31— jl m€P[l] 

where P[l] is the index set of all possible periodic orbits of length I of the graph, A m [l] e C 
is the weight of orbit number m of length I, computable from the matrix elements of U, 
and L$ [I] is the reduced action of periodic orbit number m of length /. Using this result 
together with ( |3.7|) we obtain the explicit periodic orbit formula for the spectrum in the 
form 



2 M J piL^ [l]kr, 



kn = K - -Im£y £ A^—^—sm — Lg>[i] 
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(4.10) 



Since the derivation of (|4.10|) involves only a resummation of TrS*' (which involves only 



a finite number of terms), the convergence properties of ( |4.8| ) are unaffected, and ( |4.10| ) 
converges. 

Reviewing our logic that took us from ( |4.6| ) to ( [4.10| ) it is important to stress that ( [4.1 0| ) 



converges to the correct result for k n . This is so because starting from Q4.6| ), which we 
proved to be exact, we arrive at (|4.10|) performing only allowed equivalence transformations 



(as mentioned already, the step from (|4.7f ) to (|4.8|) is proved in Sec. V). This is an important 
result. It means that even though (|4.10| ) may be a series that converges only conditionally 



(in Sec. VI we prove that this is indeed the case for at least one quantum graph), it still 
converges to the correct result, provided the series is summed exactly as specified in ( [4.101) . 



The summation scheme specified in (|4 . 1 0|) means that periodic orbits have to be summed 
according to their symbolic lengths (see, e.g., p]fl,^,|6|)p9|) and not, e.g., according to their 
action lengths. If this proviso is properly taken into account, ( f4.10|) is an explicit, convergent 



periodic orbit representation for k n that converges to the exact value of k n . 

It is possible to re-write ( |4.10| ) into the more familiar form of summation over prime 



periodic orbits and their repetitions. Any periodic orbit m of length I in ( f4.10p consists of 



an irreducible, prime periodic orbit m v of length lp which is repeated v times, such that 

l = vl v . (4.11) 

Of course v may be equal to 1 if orbit number m is already a prime periodic orbit. Let us 
now focus on the amplitude A m [l] in (|4.8|). If we denote by A mp the amplitude of the prime 
periodic orbit, then 

A m [l] = l v A» mv . (4.12) 

This is so, because the prime periodic orbit m-p is repeated v times, which by itself results 
in the amplitude A" ■ The factor lp is explained in the following way: because of the trace 
in ( |4.8| ), every vertex visited by the prime periodic orbit mp contributes an amplitude A v mv 
to the total amplitude A m [Z]. Since the prime periodic orbit is of length lp, i.e. it visits lp 
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vertices, the total contribution is l-p A v mv . Finally, if we denote by LjW the reduced action 
of the prime periodic orbit %, then 



(4.13) 



Collecting ( [4.11|) - (|4. 13| ) and inserting it into Q4.10Q yields 

2 , ^ 1 ^ 1 



k n = k n - £ Im ~W E ^ A mv e ivL ^p' kn sin 



L (0) 

2to mv 



(4.14) 



where the summation is over all prime periodic orbits m-p of the graph and all their repeti- 
tions v. It is important to note here that the summation in ( |4.14| ) still has to be performed 
according to the symbolic lengths I = vl-p of the orbits. 

In conclusion we note that our methods are generalizable to obtain any differentiable 
function f(k n ) directly and explicitly. Instead of integrating over N(k) alone in Q4.5D we 
integrate over f'(k)N(k) and obtain 



f(k n ) = nf{k n ) - {n - l)/(fcn_0 - / f'(k) N(k) dk. 

fen — 1 



(4.15) 



Following the same logic that led to ( |4.10|) , we obtain 



2 00 1 

f(k n )^nf(k n )-(n-l)f(k n ^)--lmJ2j E A m [l]G n (L^[l]), (4.16) 



1=1 I meP[l] 



where 



GJx) 



f{k)e ixk dk. 



(4.17) 



This amounts to a resummation since one can also obtain the series for kn first, and then 



form f(k n ). 



V. INTERCHANGE OF INTEGRATION AND SUMMATION 



One of the key points for the existence of the explicit formula ( [4.1 0|) is the possibility to 
interchange integration and summation according to 
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rb / 00 pina(x) \ 00 / „ft p ina{x) 

la (J! ^ = S U 



dx , (5.1) 



where cr(x) is continuous and has a continuous first derivative. We will prove (|5.1|) in two 
steps. 

Step 1: According to a well-known theorem (see, e.g., Ref. |!6|, volume I, p. 394) summa- 



tion and integration can be interchanged if the convergence of the sum is uniform. Consider 
an interval [<Ji, a 2 ] that does not contain a point ctq = mod 2tt. Define 



00 gincr 



<%) = E — ■ (5-2) 

n=l U 

Let 

m = - ^ — f — -, — rr + i 7r ~ ain0d27r . (5.3) 

M ; 2 2[l-cos(a)] 2 v ; 

Then, according to formulas F2 and F3 (Appendix B), S(a) = /(er) in [01,02]. In other 
words, S(a) is the Fourier series representation of f(cr). According to another well-known 
theorem (see, e.g., Ref. |37]], volume I, pp. 70-71) the Fourier series of a piecewise continuous 
function converges uniformly in every closed interval in which the function is continuous. 
Since f(a) is continuous and smooth in [01,02], S(a) converges uniformly in [01,02]. This 
means that summation and integration can be interchanged in any closed interval [xi,^] 
for which 0(2;) 7^ mod 2tt Va; G [x\, Xq\. 

Step 2: Now let there be a single point x* G (£1,22) with a(x*) = mod 2ir. Then, for 
any ei, e 2 > with x* — e\ > xi, x* + e 2 < x 2 , S(a(x)) is uniformly convergent in [xi,x* — e x ] 
and [x* + €2,^2] and integration and summation can be interchanged when integrating over 
these two intervals. Consequently, 

r x 2 / 00 P incr(x) \ 

L (S — ) * ■ 

rx*-ei / 00 p in<r(x)\ r%*+£2 ( 00 p in<r(x) \ fa!2 / 00 ina(x) \ 

L (S — ) + /,-„ (S — ) + U (S — ) * = 

x / /-a;*-ei p mcr(a;) \ x / -.^ jncr(a;) \ rx*+E 2 ( 00 p ina(x) \ 

£ (/, — + £ (/, +u — *°) + /,-„ (S — ) <"> 
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Since Y^=\ JTi ~ ei ^xp[ina(x)}/ndx and Y^=i Ix*+t 2 exj p[^ ncr ( x )]/ n dx are both uniformly 
convergent, we have with Lll (Appendix B): 

X2 / oo ina(x) \ 

y - — dx = 



E 

n=l 



x*-ei e ina(x) „ X2 E ina(x) 

dx + dx 



+ (E—-)dx. (5.5) 



I xi n Jx*+E2 Ti 

Since exp[ina(x)]/n is a non-singular, smooth function at x = x*, there is no problem with 
taking ei,e2 — > for the first two integrals on the right-hand side of ( |5.5| ). Therefore, 
integration and summation on the left-hand side of ( |5.5| ) can be interchanged if 

rx*+e 2 ( 00 p ina(x)\ oo / rx *+ e2 p ina(x) 



lim 



/ V dx = = lim V / dx . (5.6) 



This is guaranteed according to T3 (Appendix A). Assuming that cr(x) has only a finite 
number N of zeros mod 2ir in (a,b), we can break (a, b) into iV sub-intervals containing a 
single zero only, in which the interchange of summation and integration is allowed. This 
proves ( |5.1|) . 

Returning to the crucial step from (|4.7|) to (|4.8|) we have to show that 



rk 2 00 1 00 1 rk 2 

/ y -S n (k)dk=y - S n (k)dk. (5.7) 
Since ^(A;) is unitary, it is diagonalizable, i.e. there exists a matrix VK(fe) such that 

S(k) = W{k)di&g(e iai{k \...,e ia2N B ik) ) W\k), (5.8) 



where <Ji(k), . . . ,cr 2 N B {k) are the 2N B eigenphases of S{k). Because of the structure (p.2[) 
of the S matrix in conjunction with the smoothly varying phases ( p.3|) , the eigenphases of 
the S matrix have only a finite number of zeros mod 2tt in any finite interval of k. This is 
important for later use of ( |5.1|) which was only proved for this case. 

We now make essential use of our focus on finite quantum graphs, which entails a finite- 
dimensional S matrix, and therefore a finite-dimensional matrix W. In this case matrix 
multiplication with W leads only to finite sums. Since for finite sums integration and 
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summation is always interchangeable we have 

Jnax(k) p ina 2 N B (k) 



rk 2 00 1 rk 2 00 / „mai{k) p ™a 2 N B (k)\ r~X 

I E " * = /„ »' (*) E diag (— , .... — ^— J H-'» (*> * y 

00 fc2 / inai(k) incr 2 N B {k)\ 00 1 r fc 2 

y/ ^(^diag )W\k)dk=Y- S n {k)dk. (5.9) 

^ x Jkx \ n n J J kl 



This equation justifies the step from (|4.7Q to (|4.8| ), which proves the validity of ( [4.10|) and 

(EE)- 

VI. CONVERGENCE PROPERTIES OF THE PERIODIC ORBIT SERIES 

In this section we prove rigorously that ( |4.14| ) contains conditionally convergent as well as 
absolutely convergent cases. We accomplish this by investigating the convergence properties 
of ( gap in the case of the dressed three-vertex linear graph shown in Fig. 5. The potential 



on the bond Byi is zero; the potential on the bond B23 is a scaling potential explicitly given 
by 

U23 = \E, (6.1) 

where E is the energy of the quantum particle and A is a real constant with < A < 1. 
The quantum graph shown in Fig. 5 was studied in detail in |3-12.29]. Denoting by a the 



geometric length of the bond Bi 2 and by b the geometric length of the bond B 2 z, its spectral 



equation is given by [[8-12,29 



sin(u) k) = r sin(c t ; 1 /c), (6.2) 

where 

uj = a + f3b, iv 1 = a — (3b, r = j — ^ , (3 = y/l - A. (6.3) 

With 

7o = l/2, ai = r, 71 = 1/2, (6.4) 
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the spectral equation (|6.2|) is precisely of the form (|2.14|) . Since according to (|6.3|) a\ = r < 1, 



the regularity condition (|3.1| ) is fulfilled and ( |6.2| ) is the spectral equation of a regular 
quantum graph. This means that we can apply ( [4.14| ) for the computation of the solutions 



of ( |6.2| ). In order to do so, we need a scheme for enumerating the periodic orbits of the three- 
vertex linear graph. It was shown in PSj that a one-to-one correspondence exists between 



the periodic orbits of the three-vertex linear graph and the set of binary Polya necklaces 
p9| , |38|j39| . A binary necklace is a string of two symbols arranged in a circle such that two 
necklaces are different if (a) they are of different lengths or (b) they are of the same length 
but cannot be made to coincide even after applying cyclic shifts of the symbols of one of the 
necklaces. For the graph of Fig. 5 it is convenient to introduce the two symbols £ and 71, 
which can be interpreted physically as the reflection of a graph particle from the left (Vi) or 
the right (V3) dead-end vertices, respectively. Since strings of symbols are frequently referred 
to as words, we adopt the symbol w to denote a particular necklace. For a given necklace w 
it is convenient to define the following functions |[29|| : n-ji(w), which counts the number of TZs 
in w, nc{w), which counts the number of Cs in w, the pair function a(w), which counts all 
occurrences of 7£-pairs or £-pairs in w and the function /3(w), which counts all occurrences 
of C7Z or 7ZC symbol combinations in w. We also define the function £{w) = n c {w) + nn(w), 
which returns the total binary length of the word w, and the phase function x{ w )i defined 
as the sum of £(w ) and the number of 7£-pairs in w. We note the identity 

a(w) + (3{w) = £(w). (6.5) 

In evaluating the functions defined above, we have to be very careful to take note of the cyclic 
nature of binary necklaces. Therefore, for example, a(1Z) = 1, f3(C1Z) = 2, a(CCIZIZC) = 3 
and f3(CC1Z1ZC) = 2, which also checks (|6.5|). In addition we define the set W{1) of all 
binary necklaces of length I. 

Let us look at W{2). This set contains three necklaces, ££, C1Z = TIC (cyclic rotation 
of symbols) and TZTl. The necklace CC is not a primitive necklace, since it consists of a 
repetition of the primitive symbol C The same holds for the necklace 7171, which is a 



repetition of the primitive symbol 1Z. The necklace CTZ is primitive, since it cannot be 
written as a repetition of a shorter string of symbols. This motivates the definition of the 
set W-p of all primitive binary necklaces and the set Wp(l) containing all primitive binary 
necklaces of length I. 

An important question arises: How many primitive necklaces N-p(l) are there in W{1)1 
In other words, how many members are there in Wp(l) C W(l)7 The following formula gives 



the answer 39 



Nv{l) = j E 2' /m , (6-6) 

m\l 

where the symbol u m\V denotes "m is a divisor of Z", and 4>(m) is Euler's totient function 
defined as the number of positive integers smaller than m and relatively prime to m with 
0(1) = 1 as a useful convention. It is given explicitly by |4C 



0(1) = 1, 0(n) = n J] (l - -) , n > 2, 

pin V PJ 



(6.7) 



p\r 

where p is a prime number. Thus the first four totients are given by 0(1) = 1, 0(2) = 1, 
0(3) = 2 and 0(4) = 2. 

A special case of ( |6.6|) is the case in which I is a prime number. In this case we have 
explicitly 

N v (p) = - (2 P - 2) , p prime. (6.8) 
P 

This is immediately obvious from the following combinatorial argument. By virtue of p 
being a prime number a necklace of length p cannot contain an integer repetition of shorter 
substrings, except for strings of length 1 or length p. Length p is trivial. It corresponds 
to the word itself. Length 1 leads to the two cases 1Z1Z1Z1Z1Z... 1Z and CCCCC..C, where 
the symbols 1Z and £, respectively, are repeated p times. So, except for these two special 
necklaces, any necklace of prime length p is automatically primitive. Thus there are 



1 P 



P 



(6.9) 



v 
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different necklaces with v symbols £ and p — v symbols TZ, where the factor 1/p takes care 
of avoiding double counting of cyclically equivalent necklaces. In total, therefore, we have 

1 



N v {p) 



v 



(2 P - 2) 



(6.10) 



primitive necklaces of length p, in agreement with fl6.8|) . The sum in ( |6 . 1 0| ) ranges from 1 
to p — 1 since v = would correspond to the composite, non-primitive necklace 1Z1Z1Z1Z...1Z 
and v = p would correspond to the composite, non-primitive necklace £££££...£. 

We are now ready to apply ([4.1 4|) to the three- vertex linear graph. In "necklace notation" 
it is given by [§-[l2],|29| 

A" . 



h 



2 oo oo 

r X!X X! ,, , in] 

" 1=1 v=\ wew P -Mwew{i) v 



sin 



vL$ k n 



sm 



VTX 



L (o) 



-2^o 

where Lffl is the reduced action of the primitive necklace w, given by |29| 



(6.11) 



= 2[n c {w)a + n n {w)(3b] 



and the amplitude A w of the primitive necklace w is given by [E9| 



(6.12) 



A w = (-l)xW r a ^ (1 - r 2 )^' 2 



(6.13) 



where r and (3 are defined in ( |6.3D . The notation vw refers to a necklace of binary length 
vt{w) that consists of v concatenated substrings w. Note that the summations in ( |6.11|) are 
ordered in such a way that for fixed I we sum over all possible primitive words w and their 
repetitions v such that the total length of the resulting binary necklace amounts to /, and 
only then do we sum over the binary length / of the necklaces. This summation scheme, 
explicitly specified in ( |6.11|) , complies completely with the summation scheme defined in 
Sec. IV. Since we proved in Sec. IV that ( |4. 14| ) converges, provided we adhere to the correct 
summation scheme, so does ( |6.11| ). 



A numerical example of the computation of the spectrum of ( |6.2| ) via ( |6.11| ) was presented 
in || where we chose a = 0.3, b = 0.7 and A = 1/2. For n = 1, 10, 100 we computed the exact 
roots of (|6.2j ) numerically by using a simple numerical root-finding algorithm. We obtained 
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fc (exact) _ 4 . 10 7149 ? fci e xact) « 39.305209 and ^ e o act) « 394.964713. Next we computed 
these roots using the explicit formula (|6.11|) . Including all binary necklaces up to I = 20, 
which amounts to including a total of approximately 10 5 primitive periodic necklaces, we 
obtained « 4.105130, « 39.305212 and ifejg^ « 394.964555. Given the fact 

that in Sec. IV we proved exactness and convergence of ( |4.14j ) ( (|6.11|) , respectively), the 
good agreement between k^ xact ^ and fc^ p '° - ) , n = 1,10,100, is not surprising. Nevertheless 
we found it important to present this simple example here, since it illustrates the abstract 
procedures and results obtained in Sec. IV, checks our algebra and instills confidence in our 
methods. 

We now investigate the convergence properties of fl6.11|) for two special cases of dressed 
linear three-vertex quantum graphs (see Fig. 5) defined by 

r = — y=, a = m/3b, m—1,2. (6.14) 



In this case the reduced actions (|6.12|) reduce to 



and co>o is given by 



= 2a[n c (w) + n n (w)/m} (6.15) 



w = o(l + — ). (6.16) 

m 



Using Q6.51 ), the amplitudes Q6.13|) are 



^ = (-1)^2-^/2. (6.17) 



We now show that for m — 1 the first sin-term in (|6.11|) is always zero, and thus Q6.ll ) 
converges (trivially) absolutely in this case. For m = 1 (|6.15| ) becomes 

L$ = 2a[n c (w) + n n {w)) = 2a£(w). (6.18) 



Also, according to fl3T3] ) and (|6.4f) k n is given by 



k n = — [n + fi + 1] . (6.19) 
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Thus, for m = 1 the argument of the first sin-term in (|6.11|) is given by 



vL^kn = u£(w)(n + fx + 1 



7T. 



(6.20) 



This is an integer multiple of tt, and thus all terms in the periodic-orbit sum of (|6.11 ) vanish 
identically. Therefore we proved that there exists at least one case in which the periodic-orbit 
sum in ( |6.11 ) is (trivially) absolutely convergent. 

We now prove rigorously that there exists at least one non-trivial case in which (|6.11 



converges only conditionally. Since we already proved in Sec. IV that ( |6.11|) always converges, 
all we have to prove is that there exists a case in which the sum of the absolute values of 
the terms in fl6.ll ) diverges. In order to accomplish this, let us focus on the case m = 2 and 
estimate the sum 

1 



* = EE E 

1=1 u=l w£W v :vw&W(l) 



sin 



sin 



(6.21) 



u 2 L <S}) 2 t(v>)/2 

We now restrict the summation over all integers / to the summation over prime numbers p 
only. Moreover, we discard all non-primitive necklaces of length p, which is equivalent to 
keeping terms with v = 1 only. Observing that trivially £(w) = p for all necklaces in Wp(p), 
we obtain: 



*>E E 

p wew r (p) 



Ll 0) 2P/2 



sin 



sin 



7T 
2Un 



L19 



(6.22) 



where the sum is over all prime numbers p. For m = 2 the reduced actions are given by 



L$ =a[2n c (w)+n n (w)} 



(6.23) 



and 



3a T 2n . 
= tt, k n = — {n + n+Vj 
2 6a 



(6.24) 



We use these relations to evaluate the arguments of the two sin-functions in ( |6.22|) . We 
obtain 



L,,, ^ k r , 



n n {w)} (n + // + 1) 



(6.25) 
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and 



7T 



(6.26) 

respectively. We see immediately that all terms in ( |6.22| ) are zero if n + fi + 1 is divisible 
by 3. This provides additional examples of (trivially) absolutely convergent cases of fl6.11| ). 
In case n + \x + 1 is not divisible by 3, only those terms contribute to (|6.22 ) for which 



2nc(w) + nji(w) is not divisible by 3. Following the reasoning that led to (|6.10| ), nc(w) 
ranges from 1 to p — 1 for w G W-p(p). Then, 2nc{w) + nn(w) ranges from p + 1 to 2p — 1 
in steps of 1. Since p + 3j is never divisible by 3 for p prime and j G N, the number of 
primitive necklaces w of length p with the property that 2nc{w) + nn{w) is not divisible by 
3 is at least 



1 

V 



1 

3p 



2 p + 2 cos 



pir 



(6.27) 



K 

where the sum over the binomial coefficients was evaluated with the help of formula 0.1521 
(TJ. Therefore, with flOSQ , QQ5D , ( gg§ , QQ7D , | sin(2jvr/3)| = v^/2 for all j G N and 



m 



2ri£(w) + n-ji{w) < 2p — 1 for w G W^>(p), we obtain 

' ^ ' r 2 p + 2cos 



s > — V — 
~ 4a Y P(2p - 1) 2p/ 2 



- 3 



(6.28) 



which obviously diverges exponentially. The physical reason is that the quantum amplitudes, 
which contribute the factor 2~ P//2 in ( |6.28j ) are not able to counteract the proliferation 2 P of 
primitive periodic orbits (primitive binary necklaces) in Q6.28 ). Analogous results can easily 
be obtained for graphs with m > 2 in ( |6.14j ). 

In summary we established in this section that the convergence properties of ( |4.14| ) 
depend on the details of the quantum graph under investigation. We proved rigorously that 
both conditionally convergent and absolutely convergent cases can be found. We emphasize 
that the degree of convergence does not change the fact, proved in Sec. IV, that ( |4.14j ) always 
converges, and always converges to the exact spectral eigenvalues. 
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VII. LAGRANGE'S INVERSION FORMULA 



The periodic orbit expansions presented in Sec. IV are not the only way to obtain the 
spectrum of regular quantum graphs explicitly. Lagrange's inversion formula |4"2"] offers an 
alternative route. Given an implicit equation of the form 

x = a + wip(x), (7.1) 



Lagrange's inversion formula determines a root x* of ( |7.1|) according to the explicit series 
expansion 

x* = a+Y] !L.^-—^( X ) , (7.2) 
provided (p(x) is analytic in an open interval / containing x* and 

ViG/. (7.3) 



\w\ < 



x — a 



ip(x) 



Since ( |3.2|) is of the form (|7.1|) , and the regularity condition ( |3.1|) ensures that (|7.3| ) is 
satisfied, we can use Lagrange's inversion formula ( \I.2\) to compute explicit solutions of 

In order to illustrate Lagrange's inversion formula we will now apply it to the solution 
of (|6.2|). Defining x = w^k, the nth root of (|6.2|) satisfies the implicit equation 



x n = irn + (—1)™ arcsin[r sm(px n )}, (7.4) 

where p = uj\/ujq and \p\ < 1. For the same parameter values as specified in || and already 
used above in Sec. VI we obtain xf xa " ct ^ = u^k^™ ^ « 3.265080, x^q^ = uokf™^ w 
31.246649 and £100 = wo&ioo ps 313.986973. We now re-compute these values using the 
first two terms in the expansion ( |7.2|) . For our example they are given by 

=7m + arcsinfr sin(pTm)] <j (-l) n + r P c ^{p^ n ) I _ 

1 — r 2 sin 2 (p7rn) J 



We obtain xf ] = 3.265021 . . ., 4o = 31.246508 ... and x% = 313.986819 . . ., in very good 

, .,-1 (exact) (exact) i (exact) 

agreement with x\ , x 10 and x^ 00 . 

24 



Although both, (|4.14|) and ( |7.2|) are exact, and, judging from our example, ( |7.2|) appears 
to converge very quickly, the main difference between ( [4.14[ ) and ( |7.2| ) is that no physi- 
cal insight can be obtained from ( |7.2j ), whereas ( |4 . 1 4| ) is tightly connected to the classical 
mechanics of the graph system providing, in the spirit of Feynman's path integrals, an in- 
tuitively clear picture of the physical processes in terms of a superposition of amplitudes 
associated with classical periodic orbits. 



VIII. DISCUSSION, SUMMARY AND CONCLUSION 

There are only very few exact results in quantum chaos theory. In particular not much 
is known about the convergence properties of periodic orbit expansions. Since quantum 
graphs are an important model for quantum chaos |[43|| , which in fact have already been 
called "paradigms of quantum chaos" p5] , it seems natural that they provide a logical 
starting point for the mathematical investigation of quantum chaos. The regular quantum 
graphs defined in this paper are important because they provide the first example of an 
explicitly solvable quantum chaotic system. Moreover regular quantum graphs allow us to 
prove two important results: (a) Not all periodic orbit expansions diverge. There exist 
nontrivial, convergent, periodic orbit expansions, (b) There exist explicit periodic orbit 
expansions that converge to the exact values of individual spectral points. 

The main result of this paper is an analytical proof of the validity and the convergence 
of the explicit spectral formulas ( |4. 10| ) and (|4.14|) , respectively. This result is novel in two 
respects, (i) While periodic orbit expansions of the spectral density and the spectral staircase 
of a quantum system are basic tools of quantum chaos, the very concept of a periodic orbit 
expansion for individual spectral eigenvalues is new. (ii) Due to the exponential proliferation 
of the number of periodic orbits with their (action) lengths, it is frequently assumed in the 
quantum chaos community that periodic orbit expansions are formal tools at best, but do 
not converge. We proved in this paper that, at least as far as regular quantum graphs are 
concerned, and despite the exponential proliferation of periodic orbits in this case ||, the 
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periodic orbit expansion (|4.14j) converges in the usual sense of elementary analysis. This 
result is also new. 



The main ingredient in the proof of ( |4.14j ) is theorem T2 (Appendix A), i.e. an analytical 



proof that there is exactly one spectral point in every root cell I n . In discussions with our 
colleagues we found that while many pointed out the necessity of justifying the interchange 
of integration and summation in (|5.1|) (now established in Sec. V with T3 (Appendix A)), 
many were initially puzzled by the existence of root intervals and the organization of the 
spectral points into root cells, now guaranteed by T2 (Appendix A). This is so because 
regular quantum graphs have a positive topological entropy |2]|],f| and are in this sense 
quantum chaotic systems. Hence the spectrum of regular quantum graphs is expected to be 
"wild", in complete contrast to the fact, proved in this paper, that the spectrum of regular 
quantum graphs can actually be organized into regular root cells. In this sense regular 
quantum graphs are closely related to other quantum chaotic systems that also show marked 
deviations from the expected universal behavior [|6|,[7|,f4"4"f . As a specific example we mention 
chaotic billiards on the hyperbolic plane generated by arithmetic groups We hope that 
the pedagogical presentation of the proofs in Appendices A and B, with their hierarchical 
structure and the use of only elementary analysis concepts will help to establish theorems T2 
and T3 (Appendix A), and their consequence, the existence of explicit, convergent periodic 
orbit expansions. We mention that the spectral equation (|2.14|) of a finite quantum graph is 



an example of an almost periodic function ||46|| . More information on the analytical structure 



of the zeros of almost periodic functions can be found in [47 



There are many basic quantum mechanical problems that lead to transcendental equa- 



tions of the type (|2.14| ). So far the recommended method is to solve them graphically or 



numerically (see, e.g., [[4^,[4j|). Based on the results presented in this paper, a third method 
is now available for presentation in text books on quantum mechanics: explicit analytical 
solutions. When the regularity condition ( |3.1| ) is satisfied, either the Lagrangian inversion 
method or the periodic orbit expansion ( f4.14| ) may be employed. Since the Lagrangian 



method is a purely mathematical tool without immediate physical meaning, the periodic 
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orbit expansion may be preferred due to its direct physical relevance in terms of concrete 
classical physics. 

Having been established with mathematical rigor in this paper, formula ( 4.14| ) may serve 



as the starting point for many further investigations. We mention one: Since according 
to (4.14]) k n is known explicitly, so is the level spacing s n = k n — k n -\. This may give us 



an important handle on investigating analytically and exactly the nearest-neighbor spacing 
statistics P(s) |6|j4^| of regular quantum graphs. Whatever the precise properties of P(s) 



will be, one result is clear already: due to the existence of the root-free zones F n , established 
in Sec. Ill, P(s) is not Wignerian. Thus, regular quantum graphs will join the growing class 
of classically chaotic systems which do not show the generic properties of typical quantum 
chaotic systems. 

A corollary of some significance is the following. Since we proved that for regular quantum 
graphs there is exactly one root k n of (|2.14|) in I n , this proves rigorously that for regular 



quantum graphs the number of roots of ( |2.14| ) smaller than k grows like N(k) ~ u^k/it 
(Weyl's law §). 

An open question is the generalization of our results to the case of infinite quantum 
graphs. In case J2iZi \ a i\ converges, it seems straightforward to generalize the regularity 
condition (|3.1|) to the case of infinite quantum graphs. 

In summary we proved a rigorous theorem on the existence and convergence of explicit 
periodic orbit expansions of the spectral points of regular quantum graphs. We hope that 
this paper will lay the foundation for further rigorous research in quantum graph theory. 
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X. APPENDIX A: THEOREMS 



Theorem Tl: Let a u u h a<j G R, % G I := {1, . . . , N}, N G N, Eili |oj| < 1, and < 1. 
Then: 

2 



fix) 



E i= i Oi^j sin^x + on 



1 



Ei=l Ot COs(UiX + Qti) 



< 1 Vx G R. 



(10.1) 



Proof: Define 0j := ajjX + c^, i € I. Then: 



f(x)< 



E2=i |o»| Mil I sin(^i) | 



1 



Ei=i I cos(6»j)| 



< 



Eili H I sin(6»i)| 



Ei=i |ai| I cos(6»i)| 



(10.2) 



Define the three functions: 

N 



N 



S 2 (x) 



S(x) :=J2\ a i\ sin ( x i)^ C(x) :=^2\ai\cos(xi), g(x) := — — 

i=\ i=i 1 ~ ° l x 



(10.3) 



where x := (xi, . . . ,xji) G H N . Since there is always an x such that |sin(#j)| = sin(xj 
| cos(8i)\ = cos(xj), we prove Tl by showing that g(x) < 1 Vx G R w . Because of C 2 (x) 



\<k\ |%| cos(xj) cos(xj) < E« |%| = (E? |oi|) 2 < 1 Vx G R^, the function g is well- 



defined and singularity-free in R . Since g is differentiable in R we prove g < 1 by looking 
for the extrema of g: 

dg(x) 



dxi 



(1 - C\cc))S(a!) cos(x k ) - S\x)C(x) sin(x fe ) = 0, k E I. (10.4) 



Let x* be a solution of (|10.4|) . There are three different cases, (i) S(x*) = 0. In this case we 
have g(x*) = < 1. (ii) C(x*) = and S(x*) ^ 0. In this case we have g(x*) = S 2 (x*) < 1. 
(iii) C(x*) ^ and S(x*) ^ 0. In this case ( [101 ) reduces to 

(1-C 2 (x*)) 



S(x*)C(x*) 
For g evaluated at x* of ( |10.5| ) we obtain: 



(Os(x* k ), k =G /. 



(10.5) 



1 -C 2 (x*) 
S(x*) C(x*) 



Ej=i | Oil cos(x-) 



1 - C 2 (x*) 



1 - C 2 (x*) 
S 2 {x*) 



g(x* 



(10.6) 
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This implies g 2 (x*) = 1, or, since g > in H N , g(x*) = 1. Since there are no boundaries to 
consider where absolute maxima of g might be located, the local extrema of g encompass 
all the maxima of g in H N and we have g < 1 in R . This proves Tl. 
Theorem T2: Consider the spectral equation 

F(x) := cos(x) - $0) = 0, (10.7) 

where 

N 

$(x) = <2j cos(cG>j£ + a>i) (10.8) 

i=l 

with a,, Ui, on, x G R, % G / := {1, . . . , iV}, iV G N, \ a i\ < 1; an d < 1- Then there 
is exactly one zero x* n of ( |10.7|) in every open interval I n = {y n , v n+ i), v n = mi, n G Z. 



Proof: 

(i) First we observe that |$(x)| < J2iLi < 1 Vx G R. 



(ii) We use (i) to verify that the points v n are not roots of (|10.7|): \F{v n )\ = |(— l) n — $(z/ n )| > 



(i) 

1 — |$(z/ n )| >0. This means that roots of ( p.0.71 ) are indeed found only in the open intervals 

4- 

(hi) We define the closures I n = [u n , v n +i\- 111 In we define £ according to 

X = U n + £, < £ < 7T. (10.9) 

Inserting ( |10.9| ) into ( |10.7|) we see that in I n the spectral function F(x) is identical with 

/„(o = (-iroos(o-MO, (io.io) 

where 

N 

PniQ = ^2 a i cos(uJi£, + a { + niTUi). (10.11) 

i=l 

(iv) Because of (i): signF(i/ n ) = (—1)™. We use this fact to show: sign F{v n )F{v n +i) = 
(— l) 2n+1 = —1. Since F is continuous, this proves that there is at least one root of F in 
every I n , n G Z. 
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(v) According to (iii) and ( |10.10|) the roots of F in I n satisfy (— l) n cos(£) = <p n (Q, or 

c = a»(0, ( 10 - 12 ) 

where — arccos[(— l) n <p n (£)]. Therefore, roots of F are fixed points of (3 n . 

(vi) In I n : 

aiUi sm(ujif + a,; + nirUi) T i 

mor= [ ; — ! — %<i i. ads) 

1 — X);=i a i cos(o;j^ + «j + nnuji) 

(vii) Because of (vi) the conditions for L16 are fulfilled and f3 n has at most one fixed point 
in I n . This means that F has at most one root in I n . Since according to (iv) there is at 
least one root of F in I n , it follows that F has exactly one root in I n . 

Theorem T3: 

oo ina{x) \ 00/ „ x *^ e2 pincr(x) 



lim 

ei,e2— »0 J x * 



/ V- Ux = = lim V / dar , (10.14) 



where 61,62 > 0, a(x*) mod 27r = and a(x) is continuous and has a continuous first 

derivative. 

Proof: 

The two limits in ( |10.14j ) are independent. Therefore, splitting the integration range into 
two pieces (allowed with Lll), one from x* — ex to x*, and the other from x* to x* + €2, it 
is enough to prove 

rx*+e ( 00 p ina{x)\ oo / rx *+e p ina(x) \ 

where e > 0. The case e < 0, covering the other integral in ( |10.14|) is treated in complete 



lim 

e->0 



analogy. The first equality in ( 10.15|) follows immediately from F2, F3 and the fact that 



according to L20 the real part has a Riemann-integrable log singularity at x = x* and the 
imaginary part has a Riemann-integrable jump- singularity at x = x*. The second equality 
is more difficult to prove. 

For the following considerations we assume a'(x*) 7^ 0. We will comment on the case 
a'(x*) = below. According to L10 3e* > : f(x) = a'{x)/a'{x*) > 1/2, for \x — x*\ < e*. 
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Let e < e*: 



rx 
Jx* 



x * +e exp[ina(x)] 



dx 



n 



L10 

< 



x *+^ exp[ina{x)} g a\x) ^ 



n 2 \a'(x*) 



\exp[ina(x* + e)] — 1| < 



n cr'(x*) 
4 



n 2 |o"'(a;*)| ' 



10.16) 



where the last estimate is a simple consequence of the fact that the exponential function is 
unimodular. While this simple estimate will be useful later on for the case of large n, we 
need a better estimate for small n: 

f . \a(x* + e) - a(x*)] ) 



\exp[ina(x* + e)] — 1| 



exp < me- 



(10.17) 



Now because a(x) is different iable, we have according to the intermediate value theorem of 
differential calculus: 3£(e) : [a(x* + e) — a(x*)]/e = cr'(£); £ G [x*,a;* + e]. Therefore: 



L6_ 



|exp[mcr(x* + e)] - 1| = |exp[mcr'(£)e] - 1| <2|n<7'(£)e| = 2ne\a'(^)\ 



for |n(T'(O e l < 2, or, ne|<r'(0l/2 < 1- Let 



(10.18) 



iV(e) 



where |_ J is the floor function ([x\: largest integer smaller than x). Then: 



10.19) 



rx 
Jx* 



x * +t exp[ino"(a;)] 



dx 



n 



< 



4e 



n 



a'(x*) 



L10 6C C / AT/ \ 

< — for n < N(e). 
n 



(10.20) 



So now we have 



JV(6) 

E 

n=l 



rx 
Jx* 



n=l uX 

x * +e exp[ina(x)] 



j x * +e ^Min(r(x)} d 



n 



< 



dx 



n 



+ E 

n=N(e)+l 



rx 
Jx* 



'+ e exp[ina(x)] 



dx 



n 



< 



N(e) 



6e 



E^+ E 



„=i " „=J^ )+ i n V(< 



(10.21) 



Both sums vanish in the limit of e — > 0. We show this in the following way. For the first 
sum we obtain: 



N{e) 6e L9 / 

]T — < 6e[l + ln(7V(e))] = 6e 1 + In 

n=l n V 



1 



eK£(e))l 



< 
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6e \l - ln(e)] - In 



In 




-ln|a'(x*)|| 




a'(x*) 





< 



6e{l - ln(e) + ln(2) - In \a'{x*)\} ^ for e -> 0. 



(10.22) 



For the second sum we obtain: 
4 



E 



L8 

< 



< 



12e 



n= j% +1 n 2 |a'(x*)| ^'(x*)!^) 4-3eK£(e))| 
where we used 

2 



for e -> 0, 



10.23) 



iV(e)> 



^)) 



L10 
1 > 



3e\a'(x* 



1. 



(10.24) 



For the above proof we assumed cr'(x*) ^ 0. But our proof still works for the case a'(x*) = 
if we use a"(x)/a"(x*) instead of a'(x)/a'(x*) in ( 10.16Q . After a partial integration and 
noting that (i) lim e _ a'(x* + e) = and (ii) V|x - x*\ < e3C(e) > 0, lim e ^ Q C(e) = : 
[cr'(a;)] 2 < C(e)\a'(x)\, we arrive at an equation very similar to (|10. 16 ) . Then, following 
the steps ( |10.17|) - (|10.23| ) establishes T3 in this case too. This idea can be generalized 
to the case where the first non- vanishing derivative is of order n, i.e. a^ k '(x*) = for 
k = 0,1,. ..,ra - 1, a^(x*) ^ 0. In this case we use {x) / {x*) in (|mTB|) . It is 
not possible that all derivatives of cr(x) vanish at x* since (making the physically justified 
assumption that cr(x) is an entire function) this would mean that <j(x) is identically zero and 
there exists a scattering channel in which "nothing happens" . Since these trivial scattering 
channels have no influence on the spectrum of a given quantum graph, it is possible to 
eliminate all trivial scattering channels and thus reduce the dimensionality of the S matrix 
such that none of the eigenphases of the new, reduced S matrix is constant. Taking this 
into account, our proof establishes T3 without any exceptions. 



XI. APPENDIX B: FORMULAE, DEFINITIONS AND LEMMAS 



This appendix is a collection of formulae, definitions and lemmas needed for the proofs 
presented in the text and in Appendix A. They are collected here since they are lowest in 
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the hierarchy of proof ideas. "Formulae" are identities that can be found in tables. We 
compiled them in this appendix for completeness and easy reference. "Lemmas" are simple 
theorems of a general nature that may be found in analysis text books, but are not usually 
easily accessible. So we compiled them here for completeness and convenience. "Theorems" 
are specific to the context of this paper. They are harder to prove and unlikely to be found 
in standard analysis text books. Every lemma and theorem is proved explicitly, unless a 
convenient proof is found in the literature (see, e.g., Lll). 
Formula Fl: 1 — cos(x) = 2 sin 2 (a;/2). 



Formula F2 



Formula F3 



(see [|T1 formula 1.4411, p. 44): £~ x = ^p, < x < 2vr. 



sec 



[0 formula 1.4412, p. 44): E~i ^ = jln(5F=S5r) = ~ ln ( 2 ) 



ln[sin(x/2)], < x < 2n. 

Definition Dl: A function / is different iable in x -<=^ 3 A G R V<5 > Be* (5) > V|e| < e*: 

f( x + e )-f{x) 



A 



< S. 



The constant A is also denoted as A = f'{x). 



Lemma LI: Let a, b G R. Then: \a + ib\ = \J a 2 + b 2 < \a\ + \b\ (*). 

Proof: a 2 + b 2 < a 2 + 2\a\ \b\ + b 2 = (\a\ + |6|) 2 . Monotony of the square root yields (*). 

Lemma L2: 1 — x 2 /[n(n + 1)] > for |x| < 1 and n G N. 

Proof: (|x|/n)(|x|/(n + 1)) < 1 1 - x 2 /[n(n + 1)] > 0. 

Lemma L3: | sin(x)| < |x| Va; G R. 

Proof: Trivial for x = and \x\ > 1. Let \x\ < 1: | sin(x)| = sin(|x|) = \x\ — \x\ 3 [l — x 2 /(4 

L2 

5)]/3! - |x| 7 [l - x 2 /(8 ■ 9)]/7! - . . . < 
Lemma L4: | sin(x)| > \x\/2 for \x\ < 1. 

L2 

Proof: | sin(a;)| = sin(|x|) = \x\/2 + \x\{l - a; 2 /3)/2 + |a;| 5 [l - x 2 /{6 ■ 7)]/5! + . . . > \x\/2. 

Lemma L5: x 2 /2 < \x\ for \x\ < 2. 

Proof: x > 0: x 2 /2 - x = x(x - 2)/2 < for < x < 2. 

x < 0: x 2 /2 + x = x(x + 2)/2 < for -2 < x < 0. 

Lemma L6: | exp(ix) — 1| < 2\x\ for \x\ < 2. 
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Proof: | exp(ix) — 1| = | cos(x) +isin(x) — 1| = | — 2sin 2 (x/2) +isin(x)| < 2sin 2 (x/2) + 

L3 „ L5 

| sin(x)| < 2(a;/2) 2 + \x\ < 2\x\ for \x\ < 2. 

Lemma L7: Let f(x) be a monotonically decreasing function and f[u) = a u G R, z/ = 
1,2,.. .. Let M, G N, M < N. Then: Ef= M +i a„ < JjJ /(*) ^- 

Proof: According to the intermediate value theorem of integral calculus 3£ v G [u, v + 1] 
such that /(ar)dz = /(£„)■ Then: f£ f(x)dx = Y,»=m f" +1 f(x)dx = Y^=M /(&) > 

Lemma L8: Em+i ± < Sm £ dn = ( a -i)M^ for a > 1 - 

Proof: L7 with A" — > oo. 

Lemma L9: E„=i J < 1 + ln(iV). 

Proof: = 1 + £n= 2 £ < 1 + if > = 1 + ln(iV) . 

Lemma L10: Let f(x) be continuous in R. Then: \/x* G R with f(x*) ^ 3e*(a;*) > : 
1/2 < f(x)/f{x*) < 3/2 V|s-x*| < e*(a;*) (#). 

Proof: Since /(a;) is continuous: V<5 3e(<5) : \f(x) — f(x*)\ < 5 V|x — x*\ < e(6). Choose 
5 = 5*:= \f(x*)\/2 and define e* := e(5*). Then: | [f(x)/f(x*)\ - 1| < 5*/|/(x*)| = 
1/2 V|x — x*| < e*. This inequality can be used in two different ways: (i) 1 — f(x)/f(x*) < 
\[f{x)/f{x*)]-l\ < 1/2 f{x)/f{x*) > 1/2 V|z-a;*| < e*. This is the first inequality in 
(#). (ii) [/(*)//(**)] - 1 < \[f(x)/f(x*)] - 1| < 1/2 =► f(x)/f(x*) < 3/2 V|a; - x*\ < e*. 
This is the second inequality in (#). 

Lemma Lll: Let a„, /3„ G R, f = 1, 2, . . ., | E£li a u\ < oo, | E£li Pu\ < oo, and a, b G R. 

Define: 5 :=E~i(aa* + W- Then: |5| < oo and 5 = a (E^° =1 a,) + & (E£Li 

Proof: See ||13|| , volume II, p. 3, theorem 1.2.1. 

Lemma L12: E£li exp(iz/a;)/z/ exists and is finite in < x < 2ir. 

Proof: Follows immediately from F2, F3 and Lll. 

Lemma L13: Let f(x) be differentiable in x = a with f(a) = a and f'{a) < 1. Then 3c > a 
with c - /(c) > 0. 

Proof: Using Dl as the criterion for differentiability, we choose A = f'{a) < 1, 5 = 1 — A, 
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< 7 < €*(S) and c = a + 7. Then: 



f(a + 7) - a 



7 



- A < 



/(a + 7) -f{a) 



7 



-A 



< 5= 1 - A. 



/(a + 7) - a 



7 



< 1. 



Now: 



/(c) 



7 



/(a + 7) - a 

7 



7 > 7 



/(a + 7) - a 



7 



7 > 0. 



Lemma L14: Let f(x) be different iable in x = a with f(a) = a and /'(a) < 1. Then 3c < a 
with /(c) - c> 0. 

Proof: Analogous to the proof of LI 3. 

Lemma L15: Let f(x) be continuous in [a, b] with f(a) = a, f(b) = b, f'(a) < 1 and 
f'(b) < 1. Then there exists at least one additional fixed point z of / in [a, b] with a < z < b. 
Proof: Because of L13 there exists c > a with /(c) < c. Because of L14 there exists d < b 
with /(d) > d. Define H(x) = f(x)—x. Then H(c) < and H(d) > 0. Since / is continuous, 
so is H. Then, because of the intermediate value theorem of calculus, 3z, a < c < z < d < b, 
with H(z) = 0, i.e. f(z) = z. 

Lemma L16: Let f(x) be differentiable in [a, b] with f'(x) < 1 \fx e [a, b]. Then / has at 
most one fixed point in [a, b]. 

Proof: Assume that / has exactly n > 1 fixed points 27 < x 2 < ... < x n in [a, 6]. Then, 
because of L15, there must be at least one additional fixed point between any pair of fixed 
points Xj, j = 1, ...,n, bringing the total number of fixed points to at least 2n — 1 > n, for 
n > 1. This contradicts the assumption that / has exactly n > 1 fixed points. Therefore 
/ cannot have a finite number n > 1 of fixed points. Assume now that / has a countably 
infinite number of fixed points x±, X2, ... in [a, b]. Then, according to Weierstrafi, there exists 
an accumulation point x* in [a, b}. This implies: Ve > 3x n > x m with (x n — x m ) < e. 
Because / is differentiable we have according to Dl: 

f(x m + e) - f(x in J 



f'(Xm) 



<5 Ve<e*(5). 



(*) 
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Choose 5 = 1- f'(x m ), < e = x n - x m < e*(5), and use f(x n ) = x n , f(x m ) = x m , then 



f{x m + e) f\X m ) „// x 



f\ x n) f( x rn)_ _ ell \ 



X-n x m 
\l-f\x m )\ = l-f\x m )=5. 

This contradicts equation (*). Therefore there cannot be a countably infinite number of 
fixed points in [a, b}. Now assume that / has a continuum of fixed points in [a, b]. In this 
case we can easily show that at an interior point x* of the continuum of fixed points we have 
f'(x*) = 1 in contradiction to f'(x) < 1 \/x G [a, b]. Therefore, in summary, there cannot be 
any finite number of fixed points n > 1, nor can there be infinitely many fixed points of / 
in [a, b}. The only alternatives are zero or one fixed point, i.e. at most one fixed point, as 
stated in L16. 

Lemma L17: Let S be a unitary matrix of finite dimension B G N, B > 1. Denote by 
exp(i<7i), exp(ia 2 ), . . ., exp(icrs) its eigenvalues where <jj G R and Oj mod 2n ^ 0. Define 
M := Y.n=i S n /n. Then: |M -| < oo, i, j = 1, . . . , B. 

Proof: Since S is unitary, there exists a unitary matrix Q with S = Q diag (exp(i<7i), . . ., 
exp(ia B )) Also: S n = fidiag (exp(ino"i ),..., ex.p(ina B ))^ ■ Define the series = 
J2 < ^Li ex P(i ncr j)/ n - According to L12 these series are convergent. Then, according to Lll, 
fJ'jm '■= Hf=\ ^ji a ^^mi i s convergent, and therefore finite, because only a finite sum over a® 
is involved. Again with Lll: ^ m = £~ =1 (l/ra) E/ii fyi exp(m^)^ = E~=i(l/^(-5 n ) jm = 
Mj m . This means that since is finite, so is My, i.e. |My| < oo. 

Lemma L18: Let S be the matrix of L17. Then: for < ctj < 2ir, j = 1,...,B: 
Tr En=i S n /n = Ef=i{- ln(2) - ln[sin(a,/2)] + z(tt - ^)/2}. 
Proof: Because of Lll, 

OO 1 OO 1 OO 1 

Tr S n = E - Tr5n = E - Tr { fi diag[exp(m(7i), . . . , exp(m<7 B )] fi 1 "} = 

n=l n n=l n n=l n 

E -[exp(mcri + . . . + exp(m<7 S )] F2 ' F J ,Ln i n (2) - ln[sin((j J /2)] + i(n - <Tj)/2}. 

n=l n j=l 

Lemma L19: ImTr E^°=i S n /n = Ef=i (tt - <7,)/2. 
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Proof: Follows immediately from L18. 

Lemma L20: In | sin(x)| has an integrable singularity at x — 0. 
Proof: Let b < tt/2. Then: 

= 6 |ln(6/2) - 1| < oo. 



ro L4 ro 

/ ln[sin(x)](ix < / hn(x/2)dx 
Jo Jo 



37 



REFERENCES 

[1] R. L. Devaney, A First Course in Chaotic Dynamical Systems (Addison- Wesley Pub- 
lishing Company, Inc., Reading, Massachusetts, 1992). 

[2] E. Ott, Chaos in Dynamical Systems (Cambridge University Press, Cambridge, 1993). 

[3] R. M. May, in Dynamical Chaos, M. V. Berry, I. C. Percival, and N. 0. Weiss, editors 
(Princeton University Press, Princeton, New Jersey, 1987), p. 27. 

[4] S. M. Ulam and J. von Neumann, Bull. Am. Math. Soc. 53, 1120 (1947). 

[5] R. Bliimel and W. P. Reinhardt, Chaos in Atomic Physics (Cambridge University Press, 
Cambridge, 1997). 

[6] M. Gutzwiller, Chaos in Classical and Quantum Mechanics (Springer, New York, 1990). 

[7] H.-J. Stockmann, Quantum Chaos (Cambridge University Press, Cambridge, 1999). 

[8] R. Bliimel, Yu. Dabaghian, and R. V. Jensen, Phys. Rev. Lett. 88, 044101 (2002). 

[9] Y. Dabaghian, R. V. Jensen, and R. Bliimel, Phys. Rev. E 63, 066201 (2001). 

[10] R. Bliimel, Yu. Dabaghian, and R. V. Jensen, Phys. Rev. E, in press (2002). 

[11] Yu. Dabaghian, R. V. Jensen, and R. Bliimel, Pis'ma v ZhETF 74, 258 (2001); JETP 
Lett. 74, 235 (2001). 

[12] Yu. Dabaghian, R. V. Jensen, and R. Bliimel, JETP, in press (2002). 

[13] K. Endl and W. Luh, Analysis, three volumes (Akademische Verlagsgesellschaft, Frank- 
furt am Main, 1972). 

[14] A. S. Demidov, Generalized Functions in Mathematical Physics (Nova Science Publish- 
ers, Huntington, 2001). 

[15] W. Walter, Einfiihrung in die Theorie der Distributionen (Bibliographisches Institut, 
Mannheim, 1974). 

38 



[16] P. A. M. Dirac, The Principles of Quantum Mechanics, fourth edition (Oxford Univer- 
sity Press, Oxford, 1958). 

[17] J. von Neumann, Mathematical Foundations of Quantum Mechanics (Princeton Univer- 
sity Press, Princeton, 1955). 

[18] T. Kottos and U. Smilansky, Phys. Rev. Lett. 79, 4794 (1997); Ann. Phys. (N.Y.) 274, 
76 (1999). 

[19] S. P. Novikov, Usp. Math. Nauk, 52 (6), 177 (1997). 

[20] S. P. Novikov and I. A. Dynnykov, Usp. Math. Nauk 52(5), 175 (1997). 

[21] K. Rudenberg and C. Scherr, J. Chem. Phys. 21, 1565 (1953). 

[22] E. Akkermans, A. Comtet, J. Desbois, G. Montambaux, and C. Texier, Ann. Phys. 
(N.Y.) 284, 10 (2000). 

[23] Proceedings of the International Conference on Mesoscopic and Strongly Correlated 
Electron Systems "Chernogolovka 97", chapter 4: Quasi- ID Systems, Networks and 
Arrays, edited by V. F. Gantmakher and M. V. Feigel'man, Usp. Fiz. Nauk 168, pp. 
167-189 (1998). 

[24] H. Schanz and U. Smilansky, Phys. Rev. Lett. 84, 1427 (2000);Philos. Mag. B 80, 1999 
(2000). 

[25] U. Smilansky, J. Phys. A (Mathematical and General) 33, 2299 (2000). 

[26] F. Barra and P. Gaspard, J. Stat. Phys. 101, 283 (2000);Phys. Rev. E 65, 016205 
(2002). 

[27] P. Pakonski, K. Zyczkowski, and M. Kus, J. Phys. A (Mathematical and General) 34, 
9303 (2001). 

[28] U. Smilansky, in Mesoscopic Quantum Physics, Les Houches, Session LXI, 1994, edited 



39 



by E. Akkermans, G. Montambaux, J.-L. Pichard and J. Zinn- Justin (Elsevier Science 
Publishers, Amsterdam, 1995), pp. 373-433. 

[29] R. Blumel and Yu. Dabaghian, J. Math. Phys. 42, 5832 (2001). 

[30] Y.-C. Lai, C. Grebogi, R. Blumel, and M. Ding, Phys. Rev. A 45, 8284 (1992). 

[31] L. Sirko, P.M. Koch, and R. Blumel, Phys. Rev. Lett. 78, 2940 (1997). 

[32] Sz. Bauch, A. Bledowski, L. Sirko, P.M. Koch, and R. Blumel, Phys. Rev. E 57, 304 
(1998). 

[33] R. Blumel, P. M. Koch, and L. Sirko, Found. Phys. 31, 269 (2001). 
[34] M. Keeler and T. J. Morgan, Phys. Rev. Lett. 80, 5726 (1998). 

[35] A. C. Aitken, Determinanten und Matrizen (Bibliographisches Institut AG, Mannheim, 
1969). 

[36] R. Courant, Differential and Integral Calculus, second edition, two volumes (Interscience 
Publishers, New York, 1937). 

[37] R. Courant and D. Hilbert, Methods of Mathematical Physics, first English edition, two 
volumes (Interscience Publishers, New York, 1953). 

[38] J. H. van Lint and R. M. Wilson, A Course in Combinatorics (Cambridge University 
Press, Cambridge, 1992). 

[39] J. Riordan, An Introduction to Combinatorial Analysis (Wiley, New York, 1958). 

[40] Handbook of Mathematical Functions, edited by M. Abramowitz and I. A. Stegun (Na- 
tional Bureau of Standards, Washington D.C., 1964). 

[41] I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and Products, 6th edition, 
A. Jeffrey, Editor, D. Zwillinger, Assoc. Editor (Academic Press, San Diego, 2000). 

[42] G. Sansone and H. Gerretsen, Lectures on the Theory of Functions of a Complex Variable 

40 



(Noordhoff, Groningen, 1960). 

[43] T. Kottos and H. Schanz, Physica E 9, 523 (2001). 

[44] F. Haake, Quantum Signatures of Chaos (Springer, Berlin, 1991). 

[45] E. B. Bogomolny, B. Georgeot, M.-J. Giannoni, and C. Schmit, Phys. Rev. Lett. 69, 
1477 (1992). 

[46] H. Bohr, Almost Periodic Functions (Chelsea Publishing, New York, 1951). 
[47] B. J. Lewin, Nullstellenverteilung ganzer Funktionen (Akademie-Verlag, Berlin, 1962). 
[48] L. I. Schiff, Quantum Mechanics, second edition (McGraw-Hill, New York, 1955), p. 37. 
[49] A. Messiah, Quantenmechanik, volume I (Walter de Gruyter, Berlin, 1976), p. 89. 



41 



Figure Captions 
Fig. 1: Sketch of a quantum graph with six vertices and ten bonds. 

Fig. 2: Structure of root cell I n . To the left and to the right of I n are the root-free intervals 
and respectively. Together they form the root-free zone F n = fHufW. Roots 
of the spectral equation are found in the interval R n . The delimiters of I n are k n -\ and k n . 
The average location (star) of the root k n is given by k n . 

Fig. 3: Graphical solution of the spectral equation cos(x) = $(a;). Since cos(x) is "faster" 
than <&(x), one and only one solution exists in every interval (nn, (n + 1)tt). This fact is 
proved rigorously in Sec. III. 

Fig. 4: Detail of the spectral staircase illustrating the computation of the integral ( |4.5| ) over 
the spectral staircase from fc n _i to k n for the purpose of obtaining an explicit expression for 

Fig. 5: Three- vertex linear quantum graph. The vertices are denoted by V\, V 2 and V3, 
respectively. V\ and V3 are "dead-end" vertices. The bonds are denoted by B 12 and _B 23 , 
respectively. While there is no potential on the bond B12 (indicated by the thin line repre- 
senting the potential-free bond -B12), the bond B23 is "dressed" with the energy- dependent 
scaling potential U23 = \E (indicated by the heavy line representing the "dressed" bond 
#23) ■ 
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